Absolute dimensions of the unevolved B-type eclipsing binary 

GG Orionis^ 



Guillermo Torres 

Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138 

Claud H. Sandberg Lacy^ 
Department of Physics, University of Arkansas, Fayetteville, AR 72701 

Antonio Claret 

Instituto de Astroffsica de Andaluci'a, CSIC, Apartado 3004, E-18080 Granada, Spain 

Jeffrey A. Sabby^ 

Department of Physics, University of Arkansas, Fayetteville, AR 72701 



To appear in The Astronomical Journal, December 2000 



^Some of the observations reported here were obtained with the Multiple Mirror Telescope, a joint facility 
of the Smithsonian Institution and the University of Arizona. 

^Visiting astronomer, Kitt Peak National Observatory, National Optical Astronomy Observatories, 
operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. 



- 2 - 



ABSTRACT 

We present photometric observations in B and V as well as spectroscopic 
observations of the detached, eccentric 6.6-day double-lined eclipsing binary 
GG Ori, a member of the Orion OBI association. Absolute dimensions 
of the components, which are virtually identical, are determined to high 
accuracy (better than 1% in the masses and better than 2% in the radii) for 
the purpose of testing various aspects of theoretical modeling. We obtain 
Ma = 2.342 ± 0.016 Mq and Ra = 1-852 ± 0.025 Rq for the primary, and 
Mb = 2.338 ± 0.017 Mq and Rb = 1.830 ± 0.025 Rq for the secondary. The 
effective temperature of both stars is 9950 ± 200 K, corresponding to a spectral 
type of B9.5. GG Ori is very close to the ZAMS, and comparison with current 
stellar evolution models gives ages of 65-82 Myr or 7.7 Myr depending on 
whether the system is considered to be burning hydrogen on the main sequence 
or still in the final stages of pre-main sequence contraction. Good agreement 
is found in both scenarios for a composition close to solar. We have detected 
apsidal motion in the binary at a rate of ci; = 0.00061 ± 0.00025 deg cycle"^, 
corresponding to an apsidal period oi U — 10700 ± 4500 yr. A substantial 
fraction of this (~70%) is due to the contribution from General Relativity, 
and our measurement is entirely consistent with theory. The eccentric orbit 
of GG Ori is well explained by tidal evolution models, but both theory and 
our measurements of the rotational velocity of the components are as yet 
inconclusive as to whether the stars are synchronized with the orbital motion. 

Subject headings: binaries: eclipsing — binaries: spectroscopic — stars: 
evolution — stars: fundamental parameters — stars: individual (GG Orionis) 
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1. Introduction 

The discovery of GG Orionis (HD 290842, Tycho 4767 857 1, V = 10.4 - 11.1, 
B9.5 V, a = 05''43'"10!2, S = -00°41'15", epoch and equinox J2000) as a variable star 
is due to Hoffmeister (1934), who observed the object photographically at the Sonneberg 
Observatory. The correct period of 6.631 days was first given by Kordylewski (1951) based 
on visual and photographic times of minimum. This author obtained a mean visual light 
curve, and established that the orbit is eccentric from the displacement of the secondary 
minimum. 

Aside from the occasional measurement of the times of eclipse by a number of authors, 
GG Ori has remained until recently a rather neglected system. Double lines in the spectrum 
were detected by Lacy (1984), who described them as being narrow and of nearly equal 
strength, but no detailed spectroscopic study has been made to date. The first photoelectric 
light curves were published by Zakirov (1997), who presented light elements for this well 
detached binary solved by the method of Lavrov (1993). 

Based on the fact that the orbit is eccentric, it is expected that the system may present 
a measurable apsidal motion. This effect is of great interest in the study of detached 
eclipsing binaries because it provides information on the internal structure of stars that 
may be compared with predictions from theory. From its spectral type and other known 
properties, GG Ori was listed by Gimenez (1985) as a good candidate for the study of the 
contribution of General Relativity to the secular displacement of the line of apsides, given 
that the relativistic effect is expected to be dominant in this particular case. 

GG Ori is located in the Orion OBI association (see, e.g., Blaauw 1964; Warren & 
Hesser 1977), a complex region of star formation that has been the subject of numerous 
studies to determine the properties of the population of young stars and surrounding gas. 
The binary is located not far from the Belt of Orion, and therefore there is reason to expect 
that the system might also be quite young, adding to its interest. 

In this paper we present new high-quality photoelectric light curves in two passbands, 
which we analyze together with other published photometry. We also report the results 
of our intensive spectroscopic monitoring of GG Ori that, combined with the light curves, 
enable us to derive highly precise absolute dimensions for both components of the system. 
The stars turn out to be nearly identical in all their properties. 

Our current knowledge of the internal structure and evolution of stars is such that 
observed stellar properties determined with errors of even 5% are of little value to constrain 
the models in any significant way, since the differences between competing calculations or 
the effects of slight changes in the input physics or atomic constants are below this level 
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(see Andersen 1991, 1998). The masses we obtain here have errors smaller than 1%, and the 
radii are determined to better than 2%. We use these determinations for a comparison with 
the predictions of recent stellar evolution models, which suggest a very young age for the 
system. An analysis of all available eclipse timings along with our radial velocities leads to a 
small but apparently significant apsidal motion detection. We discuss the importance of the 
relativistic effect mentioned above compared to classical effects from rotational and tidal 
distortions. In addition, we examine the predictions from tidal theory regarding spin-orbit 
synchronization and the circularization of the orbit. 

2. Spectroscopic observations and reductions 

GG Ori was observed spectroscopically at the Kitt Peak National Observatory (KPNO) 
and at the Harvard-Smithsonian Center for Astrophysics (CfA). Observations at KPNO 
were made with the 2.1m and coude-feed spectrometer from 1983 December to 1999 March. 
A total of 15 spectrograms were obtained with a variety of CCD detectors. A spectral 
resolution of 0.02 nm or 2 pixels was used before May 1998, and 0.03 nm or 3 pixels after 
that date. The spectral coverage was about 10 nm in the first 12 spectrograms, and 32 nm 
in the last 3 spectrograms. The central wavelength of those observations was 450 nm. 

Projected rotational velocities (f sini) were determined by comparing line widths 
of the binary components' unblended features with artificially-broadened features in the 
spectra of comparison stars with known values of vsini. These comparison stars were 
chosen to be of nearly the same spectral type as GG Ori, and were observed with the same 
instrumental configuration as the binary. The reference stars used were o Peg (HR 8641, 
Al IV, vsini = 10 kms"^; Fekel 1998, private communication), 68 Tau (HR 1389, A2 IV-V, 
vsini = 18 kms^^; Hoffleit 1982), and HR 8404 (21 Peg, B9.5 Y, vsini = 4 kms"^; Fekel 
1999). From these comparisons in 10 of the spectrograms, values of 24 ± 2 kms~^ and 
23 ± 2 kms^^ were determined for the primary and secondary in GG Ori, respectively, 
based mainly on the Mg II 448.1 nm line. The uncertainties given account for the scatter of 
the line width measurements as well as the agreement between the results using different 
standard stars. The component of GG Ori we call here the "primary" (also star "A") is the 
one eclipsed at phase 0.0 in the light curve (see §4). Formally it is also the more massive 
one, but only marginally so, as we describe later. 

Radial velocities of the components were measured by cross-correlation with the 
FXCOR task in IRAFQ, using standard stars (same as above) with known radial velocities 



■^IRAF is distributed by the National Optieal Astronomy Observatories, which is operated by the 
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as templates. Radial velocities for the low-amplitude spectroscopic binary star o Peg 
were taken from the orbit of Fckcl (1999). The radial velocity of HR 8404 was taken as 
+0.2 kms^^ based on 12 velocities by Fckel (1999), and for 68 Tan we adopted the value 
+39.0 kms~^, from the same source. Our measurements for GG Ori from the 11 spectra 
with relatively unblended lines are given in Table 1. 

In addition to the radial and rotational velocities, we determined the line strength 
ratio from 17 line pairs in spectra of good signal-to-noise (S/N) ratio. This may be used 
as a proxy for the light ratio between the stars. We obtained EWa/EW^ = 1.01 ± 0.02, 
corresponding to the blue region of the spectrum. Because the two components are virtually 
identical, no correction for the difference in temperature is needed, and the ratio of the 
equivalent widths, EW^/EW^, can be assumed to be identical to the hght ratio. 

Our observations of GG Ori at the CfA were collected mostly with the cchelle 
spectrograph on the 1.5m Tillinghast reflector at the F. L. Whipple Observatory on 
Mt. Hopkins (Arizona), where the system was monitored from 1996 March to 2000 
March. Occasionally we observed with an identical spectrograph on the Multiple Mirror 
Telescope (also atop Mt. Hopkins, Arizona). A total of 42 spectra were recorded with a 
photon-counting Reticon detector covering a single echelle order centered at 518.7 nm. The 
resolution is 0.015 nm (A/AA ~ 35,000), and the spectra span about 4.5 nm. The S/N 
ratios range from about 25 to 40 per resolution element. The zero point of the velocity 
system was monitored by means of nightly exposures of the dusk and dawn sky, in the 
manner described by Latham (1992). The accuracy of the CfA velocity system, which is 
within about 0.1 kms~^ of the reference frame defined by minor planets in the solar system, 
is documented in the previous citation and also by Stefanik et al. (1999). 

Radial velocities were determined using the CfA implementation of TODCOR (Zucker 
& Mazch 1994), a two-dimensional cross-correlation algorithm that uses two templates, 
one for each component of the binary. The templates were selected from a large library 
of synthetic spectra based on the latest model atmospheres by R. L. Kurucz, computed 
for us specifically for the wavelength of our observations (Morse & Kurucz 1998, private 
communication). The instrumental profile is explicitly included by applying a Gaussian 
convolution with a FWHM of 8.5 kms^^, corresponding to our spectral resolution. These 
synthetic spectra have been calculated over a wide range of effective temperatures, projected 
rotational velocities, metallicities, and surface gravities. Initially we used templates with 
a temperature of T^q = 10000 K for both stars, surface gravity of log 51 = 4.5 (cgs), and 
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solar metallicity. The projected rotational velocities of the components were determined by 
running extensive grids of correlations against templates covering a range of values of ?;sini, 
seeking the best match to the observed spectra as determined by the correlation coefficient 
averaged over all exposures. We obtained vsini = 16 ± 1 kms^^ for both components, 
which is significantly different than our determination based on KPNO spectra. We discuss 
this further in §5 and §6. For the final templates we adopted a temperature of 9950 K 
for both stars and logg = 4.3, based on the results from the light curve solutions and 
spectroscopic orbits presented below. 

Systematic errors in the radial velocities resulting from line blending are always a 
concern, particularly when the goal is to achieve the highest possible precision and accuracy 
in the mass determinations. Although the use of a two-dimensional correlation technique 
such as TODCOR in principle minimizes those errors (see Latham et al. 1996), experience 
has shown us that this must be checked on a case-by-case basis, particularly in view of the 
narrow spectral window of the CfA observations (e.g., Torres et al. 1997; 2000). For this 
we generated artificial composite spectra by adding together synthetic spectra for the two 
components with Doppler shifts appropriate for each actual time of observation, computed 
from a preliminary orbital solution. We adopted also a light ratio close to that for the real 
stars (see below). We then processed these simulated spectra with TODCOR in the same 
manner as the real spectra, and compared the input and output velocities. The differences 
derived in this way vary systematically with phase (or radial velocity), as expected, and 
were apphed to the real velocities as corrections, even though in the case of GG Ori they 
are quite small (< 0.5 kms~^) and hardly affect the results. The final radial velocities from 
the CfA spectra with the corrections included are given in Table 2. 

The light ratio was derived from these observations using TODCOR as described by 
Zucker & Mazeh (1994). Because of the relatively small number of lines in the 4.5 nm 
spectral window observed at CfA, significant errors can be introduced in the light ratio 
when features of each component shift in and out of this region by up to 0.2 nm due to 
orbital motion. To estimate and correct for this effect we followed a procedure analogous 
to that described above for the velocities, and found the magnitude of the systematic error 
to be about 3.5%. After correcting for this we obtain La/ Lb = 1-05 ± 0.03 at a mean 
wavelength of 518.7 nm, which is sufficiently close to the visual band that we will assume it 
is the light ratio in V, given that the stars are virtually identical. 

Preliminary double-lined orbital solutions based on our radial velocities were computed 
separately for the KPNO and CfA measurements, and are compared in Table 3. Initially 
we adopted the eclipse ephcmcris given in the General Catalog of Variable Stars (GOVS) 
(Kholopov 1985): Min I (HJD) = 2,433,596.496 + 6.63147 • E. However, the residuals from 
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the orbit for the CfA data showed a systematic pattern and the scatter was significantly 
larger than expected from previous experience with similar spectroscopic material. This 
suggested that the ephemeris is not accurate enough to be propagated forward by ~ 50 yr 
(2600 cycles) to the epoch of our velocity measurements. Therefore, a new linear ephemeris 
computed from times of minimum collected from the literature was derived (see next 
section) and used here. The possibility of a systematic difference between the center-of-mass 
velocities of the primary and secondary components was explored by solving for an offset 
simultaneously with the orbital elements. This was done independently for the two data 
sets. No shift was expected because the stars are so similar. The results in the sense 
(primary minus secondary) were +1.42 ±0.97 kms~^ for KPNO and —0.74 ±0.53 kms~^ for 
CfA. Since these are indeed not significantly different from zero, the velocities for the two 
components can be considered for all practical purposes to be on the same reference frame. 

As seen in Table 3, there are shght differences in some of the elements of these two 
orbital solutions, in particular in the velocity amplitude of the secondary component (Kb), 
but we do not consider them to be significant in view of the small number of KPNO 
observations (only 11). Therefore, for the final solution discussed below we have merged the 
two data sets. 



3. Apsidal motion analysis and final spectroscopic solution 

As indicated earlier, the noticeable eccentricity of the orbit of GG Ori suggests the 
possibility of apsidal motion, even though none has been reported for this system. In 
extreme examples this can affect the spectroscopic solutions quite significantly, as in the 
case of V477 Cyg (Popper 1968). Before computing the final orbital fit, we therefore 
searched the literature for eclipse timings, both to improve the ephemeris and to investigate 
the apsidal motion. 

Table 4 lists all times of minimum available to us, of which there are 25 primary 
minima and 32 secondary minima covering about 65 yr (1930-1995). The majority were 
determined by visual or photographic means, and only the 6 recent photoelectric timings 
and one visual timing have published uncertainties. For the rest we determined the errors 
iteratively, based on the mean residuals for each type of observation from preliminary fits: 
(7pg = 0.038 days for the photographic minima, and CTy = 0.028 days for the visual minima. 
One visual estimate was found to give a large residual (0.14 days, ~ 50"), and was excluded. 
It is indicated in parentheses in Table 4. 

A linear ephemeris fit to these data gives Min I (HJD) = 2,449,717.66253(21) + 
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6.6314948(16) • E, with a phase difference between the primary and secondary minima of 
A$ = 0.42252 ± 0.00024. This is the ephemeris we used above in our prchminary orbital 
solutions, and the period is very similar to that reported by Zakirov (1997) based on a smaller 
number of observations. Although no unusual trends were seen in the timing residuals, we 
investigated the possibility of apsidal motion by computing the apparent periods separately 
for the primary and the secondary. We obtained Pj = 6.6314971 ± 0.0000021 days and 
Pii = 6.6314890 ± 0.0000030 days, which are different at the 2.2a level. Next we submitted 
these same data to an apsidal motion analysis using the method by Lacy (1992b). For 
this we assumed a fixed eccentricity equal to the average of the values derived in columns 
(2) and (3) of Table 3, and a fixed inclination angle of 89° from prehminary light curve 
solutions. The result is only marginally significant (ci; = 0.00050 ± 0.00027 deg cycle"^), 
largely due to the weak constraint provided by the early photographic and visual times of 
minimum, which have relatively low weight. Changes in those weights have little effect. 

Our radial velocities extend the time base provided by the times of minimum by 
another 5 yr, and contain valuable information on the longitude of periastron that can 
also be used. In addition, they provide a much stronger constraint on the eccentricity, 
which must otherwise be fixed when using only eclipse timings. The optimal solution, 
therefore, is to combine both kinds of measurements into a single least squares fit, solving 
simultaneously for the spectroscopic orbital elements and the apsidal motion. In this way 
the information on the period contained in the times of minimum is implicitly taken into 
account for the spectroscopic elements, rather than having to adopt a fixed ephemeris as 
we have done in §2. At the same time, the effect of the possible rotation of the line of 
apsides on the spectroscopic elements is also accounted for. In addition, we have allowed 
for an arbitrary shift between the KPNO and CfA velocity systems in view of the difference 
in the way the two zero points were established. The offset turns out to be negligible: 
(CfA-KPNO) = -0.06 ± 0.39 kms'^ 

The results of this simultaneous fit are given in column (4) of Table 3. The relative 
weights of the primary and secondary velocities in each data set, as well as the weights 
of the times of minimum, have been iterated until reaching convergence for a reduced 
of unity. Tests show once again that the weights assigned to the visual and photographic 
timings, or to the velocities, have little effect on the results. The velocity residuals from this 
solution are listed separately in Table 1 and Table 2, and the timing residuals are given in 
Table 4. The final spectroscopic orbital solution is represented graphically in Fig. 1, along 
with the observations. 

The apsidal motion resulting from the final fit is d; = 0.00061 ± 0.00025 deg cycle"^, 
which is significant at the 2.5a level, and differs from our preliminary estimate above by 
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less than half of its uncertainty. Also hsted in Table 4 are the sidereal and anomalistic 
periods (Pg and Pa), as well as the apsidal motion period U = 10700 ± 4500 yr. Pcriastron 
passage occurs at a photometric phase of 0.06. A plot of the — C deviations of the times 
of minimum from the linear terms of the apsidal motion is given in Fig. 2 along with the 
predicted deviations. The top panel draws attention to the fact that the observations 
cover less than 1% of a full cycle, while the bottom panel expands the region around 
the observations, with the coverage provided by the radial velocity measurements being 
indicated in the upper right. 

4. Photometric observations and light curve solutions 

Differential light curves were obtained by CHSL at the Cerro Tololo Inter- American 
Observatory (CTIO) with the 24-inch Lowell telescope in the Johnson B and V bands 
during 1993-1995. Absolute indices were also measured separately based on measurements 
that were carefully tied to the standard UBV system through observations of secondary 
standards from Landolt (1973), typically 30-40 standards per night. The absolute 
photometry was made at Mount Laguna Observatory near San Diego (California) in the 
fall season of 1989, and at CTIO in 1988-1990 and 1993-1995. The procedures that were 
used are described in detail by Lacy (1992a). The UBV indices on the standard system 
are shown in Table 5, and the 257 differential BV observations are given in Table 6 and 
Table 7. The differential magnitudes are in the sense (variable— comparison) and are 
referred to BD — 1°1013 (HD 38165, B9 V). The comparison stars (see Table 5) were found 
to be constant at a level of about 0.008 mag for the standard deviation of the differences 
between comparison stars. The precision of the CTIO differential observations, based on 
previous results, is estimated to be about 0.006 mag in both B and V at the magnitudes of 
the program stars. 

In addition to the absolute photometric indices mentioned above, which we have 
computed separately for the two observing intervals, uvby indices are available from Hilditch 
& Hill (1975). The results from the three different sources are somewhat inconsistent. The 
V magnitude is slightly brighter by about 2.5% in the 1988-1990 results compared to the 
1993-1995 resuhs and the vahie of Hilditch & Hill (1975), which is V = 10.380 ± 0.015. 
The B — V color index during 1993-1995 is considerably bhicr than the value measured in 
1988-1990. These discrepancies cannot be explained by differences in sky conditions (which 
were always photometric), nor by differences in aperture size, centering errors, or seasonal 
variations. Furthermore, the scatter from the light curve fits described later in this section 
is significantly larger than expected based on our previous results with the equipment used 
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at CTIO at these magnitude levels. The standard deviation of the residuals from the fitted 
photometric orbit of FS Mon (Lacy et al. 2000), for example, was 0.006 and 0.005 mag in 
B and V, respectively, whereas the corresponding figures for GG Ori (see below for details) 
are 0.014 and 0.009 mag. Both binaries were observed during the same runs. Intrinsic 
variability of one or both of the components of GG Ori is therefore suspected. 

Differential light curves in the Johnson U, B, V, and R bands have also been obtained 
in Uzbekistan by the group led by M. Zakirov. They used 0.6m telescopes at the Maidanak 
Observatory, in the southern part of Uzbekistan. Their differential magnitudes in the sense 
(variable— comparison) relative to the star BD — 1°1013 have been described by Zakirov 
(1997). The transformation of these measurements to the standard system is somewhat 
less secure than for our own observations, but there are marginal indications that the 
B — V index of GG Ori became bluer by a few hundredths of a magnitude in the interval 
1993-1994, toward the end of their observations. Though far from being conclusive evidence 
of intrinsic variability, we note that this is in the same direction as the trend suggested 
in Table 5. We have re-analyzed the differential observations by Zakirov with the same 
methods we use for our own measurements. 

The hght curve fitting was done with the NDE model as implemented in the code 
EBOP (Etzel 1981; Popper & Etzel 1981), and the ephemeris adopted is that of §3. The 
main adjustable parameters are the ratio of the central surface brightness of the secondary 
star (Jb) in units of that of the primary, the relative radius of the primary (r^) in units of 
the separation, the ratio of the radii {k = tb/ta), the inclination of the orbit (i), and the 
geometric factors e cos uj and e sin uj which account for the orbital eccentricity. As usual, we 
have allowed also for a photometric scale factor and a phase shift in all solutions. Auxiliary 
quantities needed in the analysis include the gravity-brightening coefficients, for which we 
adopt the values 1.00, 0.85, 0.70, and 0.59 in [/, B, V, and R, respectively, from Martynov 
(1973). For the refiection coefficients we adopted the value 1.0, as appropriate for stars with 
radiative envelopes. Limb-darkening coefficients (?;,) were included as one of the variables to 
be fitted, although in one case (noted below) it was fixed to an appropriate value from the 
tables of Wade and Rucinski (1985). The mass ratio [q = Mb/ Ma = 0.9982) was adopted 
from the spectroscopic analysis in §3. 

Preliminary solutions revealed a few outliers with residuals greater than ba (one of the 
CTIO measurements in B, two in V , and one of Zakirov's measurements in each of his four 
bands), which were given zero weight in subsequent iterations. 

The fitting procedure converged in the general solutions of the CTIO light curves with 
all variables adjusted simultaneously (see columns 2 and 3 in Table 8). As is often the case 
in the analysis of partially eclipsing binary stars with nearly equal components, the fitted 
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values of the parameter k, the ratio of radii, are rather uncertain, but near unity. Attempts 
at general solutions of the Zakirov light curves failed to converge. In cases such as this light 
ratios estimated from spectrograms are very important for constraining the photometric 
solutions. For GG Ori we have at our disposal two spectroscopic determinations of the light 
ratio, one approximately in the B band (from KPNO) and the other close to the V band 
(from CfA). A grid of EBOP solutions was run on the CTIO light curves for fixed values 
of k between 0.94 and 1.10, in steps of 0.01. Within this range of about 15%, the rms 
residuals from the fits were found to change very httle as a function of k (< 0.00005 mag 
in B and < 0.00006 mag in V), supporting our concerns about the indeterminacy of k 
from the light curves alone. The sum of the relative radii changes by less than 0.7% over 
the entire range of k values. The light ratios in B and V were computed from each of 
these solutions, and then by interpolation we determined the k values that correspond to 
our measured spectroscopic light ratios: k = 0.996 ± 0.011 from the KPNO spectra, and 
k = 0.972 ± 0.015 from the CfA spectra (sec Fig. 3). The value we adopt is the weighted 
average, (k) = 0.988 ± 0.009. For comparison, the average ratio of the radii derived from 
the CTIO hght curve solutions without the constraint from spectroscopy is A; = 1.03 ± 0.04 
(see Table 8), which is different by lu. 

With the ratio of the star sizes fixed to the weighted average given above, we repeated 
the light curve fits separately for the CTIO and Zakirov (1997) photometry, and the adopted 
solutions are given in Table 8 and Table 9. Attempts to include third hght as a parameter 
showed that, within uncertainties of about 1%, it was not significantly different from zero. 

The CTIO observations in B and V and the best fit models are shown in Fig. 4, and 
the data by Zakirov (1997) in UBVR along with the corresponding fitted hght curves are 
shown in Fig. 5. Table 10 gives the weighted averages of the light elements in each band 
obtained from the two data sets. In addition to and k, we list also the sum of the relative 
radii, r — rA + tb, and the radius of the secondary, r^. The errors for these quantities were 
derived as described in the Appendix. There is good agreement between the photometric 
determinations of e cos u and e sin uj and the corresponding spectroscopic values, which are 
-0.1200 ± 0.0017 and +0.1865 ± 0.0020, respectively 

The results in Table 10 show that the components of GG Ori are nearly indistinguishable 
in size and of virtually the same luminosity, which is consistent with the fact that they 
are also very similar in mass (§3). The departure from the spherical shape is insignificant. 
Both eclipses are partial, with approximately 91% of the light of the primary blocked at 
phase 0.0 (which is a transit), and 90% of the light of the secondary eclipsed at the other 
minimum. The earlier study by Zakirov (1997) reported the primary echpse to be total 
(larger secondary star blocking the primary), although the ratio of the radii was much 
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closer to 1.0. 

The uncertainties assigned to the adjusted quantities in Table 10 are generally 
the internal errors produced by EBOP. As is well known, the formal uncertainties in 
least-squares solutions in which one or more quantities have been held fixed are typically 
too small because correlations between the elements are artificially eliminated (as we 

have done by fixing k, for example). The quantities k, r^, and i are the most important 
for determining the absolute dimensions of the components. The error in k is based on 
the uncertainty in the spectroscopic light ratios, and therefore does not suffer from this 
shortcoming. For and i we have increased the formal uncertainties to account for all 
contributions to the error. 

The rms deviations given in Table 8 for the CTIO observations are somewhat larger 
than expected, based on previous experience, and suggest that one or both components of 
GG Ori may be intrinsically variable. We have searched for patterns in the residuals by 
computing the power spectra in the B and V bands, but no significant periodicities were 
found. 



5. Absolute dimensions 

The combination of the spectroscopic results in Table 3 and the light curve results in 
Table 10 leads to the absolute masses and radii for GG Ori, shown in Table 11. The masses 
are determined with a precision of 0.7%, and the radii are good to about 1.4%. 

The fact that the surface brightness ratio Jb is near unity over a wide range 
of wavelengths indicates that the effective temperatures of the stars are essentially 
the same. From Jb in the V band the difference in visual surface brightness is 
AFy = 0.0002 ± 0.0004 (Popper 1980), which translates into a completely neghgible color 
difference A(B — V) — 0.0005 ± 0.0009, and a formal temperature difference of only 8 K. 

The discrepancies in the absohite photometry at different epochs pointed out earlier in 
§4 complicate the derivation of the effective temperature of the components. Because the 
light variations may be real, we have chosen to adopted the straight average of the largest 
and smallest UBV indices for GG Ori in Table 5: V ^ 10.372±0.013, B-V = 0.511±0.021, 
U — B — 0.324 ± 0.004. The errors for the first two of these quantities are simply half of the 
difference between the 1988-1990 and 1993-1995 measurements. The reddening was derived 
using the reddening-free index Q = (U — B) —0.72{B — V) (Johnson & Morgan 1953) and the 
calibration by Deutschman, Davis & Schild (1976). We obtained {B — V)o = —0.036 ±0.006, 
{U-B)o = -0.070 ± 0.020, E{B-V) = 0.547 ± 0.022, and Ay = 1.696 ± 0.068 (using 
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Av — 3.1 • E{B — V)). The formal errors do not account for uncertainties in the cahbrations 
themselves, which are difficult to quantify. These color indices correspond to the combined 
light, but since the components are virtually identical they are also the indices of the 
individual stars. 

An independent measure of the reddening may be obtained from the uvby photometry 
by Hilditch & Hill (1975): b-y = 0.408 ± 0.015, ci = 1.002 ± 0.030, mi = +0.016 ± 0.020. 
The calibration by Crawford (1978) for B-type stars leads to {h—y)^ = —0.030 and 
E{h—y) = 0.438, from which E{B — V) = 0.592. This is slightly larger than the reddening 
derived from the UBV photometry. However, GG Ori is near the end of the range of 
validity of the uvby calibrations, where they become rather uncertain. Also, we note that 
the {b—y)o index corresponds to a spectral type of B8.5 according to Table 1 by Popper 
(1980). But the fact that we see no sign of the He I 447.1 nm hue in our KPNO spectra of 
the object strongly suggests that the spectral type cannot be earlier than B9.5. We have 
therefore chosen to rely only on the Johnson photometry above. 

The de- reddened B — V index corresponds to a temperature of 9950 ± 200 K, and a 
spectral type of B9.5 (Popper 1980). Further properties of the stars are hsted in Table 11, 
including the mean density (p), the absolute visual magnitudes, and the distance. The latter 

two are based on the visual surface brightness parameter Fy = 3.974 ± 0.003 derived from 
the intrinsic color {B — V)o and the tabulation by Popper (1980), and are thus independent 
of bolometric corrections. 

The distance modulus of GG Ori, m-M = 8.20 ± 0.10 {d = 438 ± 20 pc), is similar to 

other determinations for various subregions of the Orion OBI complex (Warren & Hesser 
1977; Brown, de Geus, & de Zeeuw 1994; Brown, Walter & Blaauw 1999; de Zeeuw et al. 
1999), and lends support to its association. 

Also listed in Table 11 are the projected rotational velocities expected if the axial 
rotations were synchronized with the mean orbital motion of the binary and with the 
orbital motion at periastron (pseudo-synchronization). The two independent measurements 
of vsini in §2 based on our KPNO and CfA spectra disagree, with the KPNO results 
being closer to the pseudo-synchronous values, while the GfA determinations suggest 
synchronization with the mean orbital motion. We discuss this further in §6 in the context 
of tidal evolution theory. 

No metallicity determination is available for GG Ori. A rough estimate can be obtained 
from our CfA spectra following a procedure analogous to that used to derive vsini, and the 
result is [m/H] = —0.15 ± 0.20, consistent with the solar abundance. 
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6. Discussion 

In this section we compare the properties of GG Ori as Usted in Table 11 with 
predictions from models. The three aspects of theory we focus on are stellar evolution, 
internal structure, and tidal evolution. For each of these we require calculations of the basic 
properties of the components, which we have taken from evolutionary tracks computed 
specifically for the exact masses we determined for the stars. The evolutionary code 
we used is that by Claret (1995). Further details on the input physics are described 
by Claret & Gimenez (1992). Convection in these models is treated with the standard 
mixing- length prescription, with a fixed mixing-length parameter of 1.52Hp that gives the 
best fit between a solar model and the observed properties of the Sun. A moderate amount 
of core overshooting is assumed (aov = 0.20i/p), although in the case of GG Ori this effect 
is insignificant due to the unevolved nature of the system. 

6.1. Stellar evolution 

Consistent with the young environment surrounding the object (the Orion OBI 
association), we find that GG Ori is indeed very close to the Zero Age Main Sequence 
(ZAMS). This is illustrated in Fig. 6, which shows a greatly enlarged section of the 
evolutionary tracks as they approach the main sequence. The full extent of the main 
sequence band is shown in the inset. The binary is so close to the ZAMS that, given the 
uncertainties, it is difficult to determine whether the stars are already burning hydrogen 
in their cores, or whether they are still in the final stage of contraction towards the main 
sequence. The error bars shown reflect the uncertainties in the effective temperatures of the 
components and their surface gravities. The error in the placement of the tracks due to the 
uncertainty in the measured masses is indicated in the lower left corner. 

The models plotted in Fig. 6 correspond to the case where the system is already on 
the main sequence (MS). Tracks for two different chemical compositions {Z, Y) are shown, 
giving metallicities very close to solar ([m/H] = +0.03, dY/dZ = 1.65, solid lines; and 
[m/H] = +0.05, dY/dZ = 1.95, dashed lines) consistent with the indication in §5. The 
enrichment laws (dY/dZ) implied by the best fit values of the helium abundance {Y) are 
within the range found in other determinations (see, e.g., Peimbert 1993; Renzini 1994; 
Pagel & Portinari 1998; Izotov & Thuan 1998), and are also similar to the values favored 
by other echpsing binaries (Ribas et al. 2000), using the same set of evolutionary models. 

Our attempts to find evolutionary tracks that fit the location of GG Ori assuming 
it is nearing the end of the pre-main sequence (PMS) phase also gave acceptable results 
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(see Fig. 7), considering the error bars, although the agreement is formally not as good as 
the main-sequence case. The PMS tracks shown in the figure are for [m/H] = +0.07 and 
dY/dZ = 2.09. 

The evolutionary ages derived for the components are quite different in the two 
scenarios (see Table 12). In the main sequence case the average age for the two stars is 
82 Myr (for Z = 0.020, Y = 0.2729) or 65 Myr (for Z = 0.021, Y = 0.2810), and the 
difference in age between the components is similar to the formal uncertainty, which is 
derived from the error in logg. On the other hand, if GG Ori is on the PMS the models 
give an average age of only 7.7 Myr. Few cases arc known of eclipsing binaries with 
well- determined absolute dimensions in which at least one component is a bona-fide PMS 
star: EK Cep (Popper 1987; Claret, Gimenez, & Martm 1995), TY CrA (Casey et al. 1998), 
and possibly RS Cha (Clausen & Nordstrom 1980; Pols et al. 1997; Mamajek, Lawson, & 
Feigelson 2000). Such systems are particularly valuable to test models of PMS evolution, 
where theory remains essentially unchallenged by observations so far. The stars in GG Ori 
are so similar and so close to the ZAMS, however, that the constraint on PMS evolution 
is not very strong. Unevolved systems such as this are much more useful for testing 
opacity and metallicity effects in the models when coupled with an accurate spectroscopic 
determination of the metal abundance, which unfortunately is not yet available for GG Ori. 
Whether or not the system has already settled on the main sequence, it is undoubtedly 
young. In this connection we note, incidentally, that light variations such as those hinted at 
in §4 are not entirely unexpected. 

It is a remarkable coincidence that among the eclipsing systems with the best known 
absolute dimensions, no fewer than six have at least one component with virtually the same 
mass as the stars in GG Ori, to within 1%. These systems are V451 Oph, YZ Cas, /3 Aur, 
WX Cep, SZ Cen (Andersen 1991), and V364 Lac (Torres et al. 1999). Thus, a total of eight 
stars (including GG Ori) can in principle be compared with a single evolutionary track, and 
they span the entire main sequence band from the ZAMS (GG Ori) to the shell hydrogen 
burning phase (SZ Cen, age ~800 Myr). In Fig. 8 wc show this comparison against one 
of the models used above for GG Ori (Z = 0.020, Y = 0.2729), for a composition close to 
solar. The agreement with theory is very good, indicating that all these systems are well 
represented by a single metallicity. 

6.2. Internal structure and General Relativity 

Our detection of the apsidal motion of the binary in §3 provides the opportunity to test 
models of the internal structure of stars. In addition, the relatively short orbital period and 
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large masses of the components along with the fact that the orbit is eccentric led Gimenez 
(1985) to propose GG Ori as a good candidate for the study of the general-rclativistic 
(GR) contribution to the apsidal motion. At the current level of the uncertainties in the 
measurement of apsidal motion in binaries, the theoretical GR contribution is separable 
from the classical (Newtonian) terms due to the gravitational quadrupolc moment induced 
by rotation and tides. The total apsidal motion can then be expressed as totot = '^n + ^gr- 
The two effects cannot be measured separately, though, so we have chosen here to compare 
the total theoretical value with the measured quantity. 

In the case of GG Ori, the GR term is predicted to be about 2.5 times larger than 
the classical term, contributing ~70% to ci;tot- Following Gimenez (1985), we obtain 
UGR = 0.000454 ± 0.000002 deg cycfe-^ 

To compute uon we have used the internal structure constants from theory (including the 
higher order terms and /C4), which are calculated at each point along the evolutionary track 
of each component. The contribution of the rotational distortions to the Newtonian apsidal 
motion rate d;jv depends on the ratio between rotational angular velocity and the Keplerian 
angular velocity. This ratio is often derived under the assumption that the stars are already 
synchronized at periastron. For a system as young as GG Ori this may not necessarily be the 
case, and in fact the observational evidence we have is conflicting (§5; see also §6.3). We have 
therefore relied on our own measurements of the projected rotational velocities, wsini, even 
though our two independent determinations from KPNO and CfA spectra are somewhat 
different (§2). The KPNO values (wAsini = 24 ± 2 kms^""^, f^sini = 23 ± 2 kms""*^) lead to 
ujn = 0.000209 ± 0.000012 deg cycle"^, while the CfA measurements {vsmi = 16± 1 kms~^ 
for both components) give ujn = 0.000183 ± 0.000009 deg cycle^^ 

When adding this to the GR term, we obtain ujtot = 0.000663 ± 0.000012 deg cycle^^ 
(KPNO) and Cotot = 0.000637 ± 0.000010 deg cycle^^ (CfA), in which the errors account for 
all contributions from measured quantities. These predictions are to be compared with the 
observed value, ujobs — 0.00061 ± 0.00025 deg cycle"^ (§3). There is good agreement with 
theory, although the uncertainty in the observed value is large enough that the constraint 
on the models is very weak. Further measurements of times of minimum are required to 
improve uJobs- 

6.3. Tidal evolution 

Tidal forces within a binary tend to synchronize the rotation of each component with 
the orbital motion, and to make the orbit circular. The two main mechanisms proposed 
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to describe these effects (Tassoul & Tassoul 1997, and references therein; Zahn 1992, and 
references therein) make somewhat different predictions for the timescales of these processes. 
In this section we compare our observations with both, ahhough the hydro dynamical 
mechanism by Tassoul has often been found to be too efficient (see, e.g.. Claret, Gimenez, 
& Cunha 1995). Synchronization times and circularization times have been computed 
as described in the latter reference and also by Claret & Cunha (1997), by integrating 
the differential equations describing the evolution of the rotation and eccentricity. All 
time-dependent properties of the stars were interpolated directly from the evolutionary 
tracks for each component. 

In performing these calculations for young stars such as GG Ori a number of 
complications arise regarding the pre-main sequence phase where the radii of the stars are 
much larger than on the main sequence. One of them is that in principle the orbital period 
(or semimajor axis) is also changing during this phase. For example, if the period of GG Ori 
is assumed to be fixed at the current value of 6.6 days, we find that the size of the orbit 
is too small to accommodate the very large radii predicted for the stars during the early 
stages of contraction along the Hayashi tracks. Thus the equations for the orbital evolution 
(period and eccentricity) and rotational evolution are coupled and must be integrated 
simultaneously (see, e.g., Duquennoy, Mayor, & Mermilliod 1992). 

Zahn & Bouchet (1989) have done this for late type stars (0.5 Mq < M < 1.25 Mq) 
and concluded that all binaries with periods up to 7 or 8 days arrive on the main sequence 
with their orbits already circularized, and that further evolution of the eccentricity on the 
main sequence is negligible. However, eccentric systems on the main sequence with periods 
shorter than this are actually observed (see, e.g., Mathieu et al. 1992). For earher-type 
binaries such as GG Ori the efficiency of the tidal mechanisms may be different, and a 
detailed investigation of the evolution in the PMS phase is beyond the scope of this paper. 

Tidal theory is still largely in development, and many aspects of it remain somewhat 
controversial. Nevertheless, useful comparisons with the observations are still possible 
under certain simplifying assumptions, and have indeed been made in a number of cases. 
Verbunt & Phinney (1995), Claret, Gimenez, & Cunha (1995), Claret & Cunha (1997), and 
others, have compared the predictions of the main mechanisms mentioned above with the 
properties of detached eclipsing systems with well-determined absolute dimensions. Their 
calculations were done in detail (i.e., by integrating the corresponding differential equations 
rather than simply using the timescales for syncronization and circularization), but avoided 
the PMS problem by excluding it altogether (Verbunt & Phinney 1995) or by including only 
the final stages of contraction, and assumed also that the orbital period does not change 
significantly, which may be a valid approximation when the stars are already close to the 
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ZAMS. 

We have followed a similar prescription here and considered only the final loop in the 
evolutionary tracks preceding the ZAMS (see Fig. 8). We start our integrations at the 
onset of the nuclear reaction ^^C(p, 7) ^^N(/3 + u) ^^C(p, 7) ^^N. At this stage this is not 
the only source of energy, since contraction is still taking place. The synchronization and 
circularization times {tsync, tare) arc not strongly dependent on the chemical composition 
in the range we have considered above. Using the formalism by Tassoul we obtain 
logtsync = 6.875 for both components of GG Ori. The predictions from the mechanism by 
Zahn, on the other hand, give logt^„^ = 8.888 and logt^„^ = 8.890, which arc significantly 
longer. The circularization times are nearly identical in both theories: tdrc = 8.888 
(Tassoul) and tdrc = 8.889 (Zahn). These values are to be compared with the mean age of 
the system counted from the starting point of the integrations, which differs by 5.4 Myr 
from the true ages given in §6.1. The modified age is then logt — 7.884 (Z — 0.020, 

Y = 0.2729) or logt = 7.775 {Z = 0.021, Y = 0.2810), assuming GG Ori is already burning 
hydrogen on the main sequence. 

We illustrate this in a shghtly different way in Fig. 9, where we focus on logg' 
as a sensitive measure of stellar evolution. The value of the surface gravity at which 
synchronization or circularization is achieved (log g^crit) is shown as a function of orbital 
period for the two main tidal theories. Stars evolve upwards in this diagram. The measured 
values of log^f for the primary and secondary of GG Ori are represented by crosses, with 
the error bars being of the same size as the symbols. According to the theory by Zahn the 
stars should not yet have reached synchronous rotation, which agrees with the measured 

V sin i values from our CfA spectra. Also, the orbit should not yet be circular, as is indeed 
observed. The theory by Tassoul also predicts that the orbit should still be eccentric, but 
indicates that the stars should already be rotating synchronously with the orbital motion, 
which appears to be more consistent with the vsini values from our KPNO spectra. 
Although as mentioned above the mechanism by Tassoul is often found to be too efficient, 
a more accurate measurement of the projected rotational velocities in GG Ori could shed 
more light on this issue. 



7. Conclusions 

New photometric and spectroscopic observations of the eccentric binary GG Ori 
combined with a reanalysis of data from the literature have allowed us to derive definitive 
orbital parameters and physical properties of the components. Our determinations have 
formal errors smaller than 1% in the masses and smaller than 2% in the radii. GG Ori 
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thus joins the elite of stars with the best estabhshed absolute dimensions. The system is 
very young, and in fact it is so close to the ZAMS that within the observational errors the 
comparison with current stellar evolution models is unable to tell us whether it has already 
settled on the main sequence (age ~ 65-82 Myr, depending on the metallicity), or whether 
it is still at the end of the contraction phase (age = 7.7 Myr). We find good agreement with 
both scenarios for a chemical composition close to solar. 

The tentative detection of apsidal motion in GG Ori with a period of U = 10700 yr has 
allowed us to make an initial test of interior structure models. This system is particularly 
interesting in that the contribution from General Relativity should be substantial (~70% of 
the total apsidal motion). Although the agreement with theory is good, further observations 
to improve the error in co are necessary for a more stringent test. 

Current mechanisms that describe the tidal evolution of binary properties disagree as 
to whether the rotation of the components of GG Ori should be synchronized with the 
orbital motion. Unfortunately our own measurements of v sin i have not settled the issue. 
Although this area of theory is perhaps the weakest of the comparisons we have discussed, 
the models do explain the significant eccentricity of the orbit system at the relatively short 
period of 6.6 days. 

We thank P. Berhnd, M. Calkins, D. W. Latham, A. Milone, and R. R Stefanik, who 
obtained many of the spectroscopic observations used here, and R. Davis, who maintains 
the CfA database of radial velocities. Wc are also grateful to Daryl Willmarth at KPNO for 
assistance with the spectroscopic observations there. Paul Etzel provided helpful comments 
on aspects of the light curve solutions. We thank the referee for useful suggestions. This 
research has made use of the SIMBAD database, operated at CDS, Strasbourg, France, and 
of NASA's Astrophysics Data System Bibliographic Services. 

8. Appendix 

The determination of the error in the relative radius of the secondary component, tb, 
in eclipsing binaries has often been a source of confusion among users of EBOP. Results 
are occasionally published in which that uncertainty is derived by propagating errors in 
the expression = k ■ ta, where k (defined as tb/ta)-! the primary radius r^, and their 
associated uncertainties are usually adopted directly from the output of the program. Aside 
from the fact that the formal errors in the latter two quantities are typically underestimated 
(see §4), this procedure is incorrect, as already pointed out clearly in the documentation 
for EBOP (Etzel 1980). The reason is that k and ta are, as a rule, strongly (negatively) 
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correlated, and simply adding their uncertainties by quadratures does not account for this 
correlation. It is still possible to make correct use of the expression above if the correlation 
is explicitly taken into account by adding the corresponding term from the off-diagonal 
elements of the covariance matrix in the least squares solution. However, this information 
is not readily available in the standard output from EBOP. Alternative procedures for 
estimating arc discussed in the documentation for the code (see also Popper & Etzel 
1981). For example, when k is fixed in the least squares solution one may use as an 
approximation the expression 

arg^k-ar^, (1) 

although this explicitly ignores the error in k, which in our case is well known independently 
of EBOP. 

The reason for the strong correlation between k and ta is that the sum of the radii, 
r = rA + tb, which is directly related to the form of the light curve, is usually very well 
determined, and therefore if k is increased the primary radius ta must decrease to maintain 
the sum. This suggests that a more sensible way to approach the problem might be to 
consider the expressions 

in which the independent variables are r and k, which are not strongly correlated in 
general. For GG Ori wc have shown this to be the case in §4. Propagation of errors in these 
equations may therefore be expected to give more reasonable estimates of cr^^ and (7^^: 
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Note that will be smaller or larger than depending only on the value of A;, as in 
eq.(l). In fact, setting 0"^ = and combining the expressions above leads back to (1). 

Practical application of these equations requires an estimate of 0"^. Although r is not 
one of the adjustable variables in EBOP, it is not difficult to obtain an estimate of its error 
by experimenting with the light curve solutions and assessing the sensitivity of r to the 
various input quantities, as we have in fact done with ta- However, to ensure consistency 
with the error in r^, we have chosen to use the first of the equations in (3) to solve for (7^, 
since cr^^ is already known and so is (7^. Once we determined the value of needed to 
reproduce 0"^^ , we used it in the expression for ■ 
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Fig. 1. — Radial- velocity curves of GG Ori from our final spectroscopic solution that includes 
times of minimum in the least-squares fit. The filled circles are for the CfA observations of 
the primary, the open circles are for the secondary, and the asterisks represent the velocity 
measurements from KPNO. 
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Fig. 2. — a) Ephemeris curve for GG Ori along with all available times of eclipse (solid 
line for the primary); b) enlargement to show the distribution of the measurements. The 
time distribution of the spectroscopic observations included in the apsidal motion analysis 
is indicated by the horizontal line in the upper right. 



-27- 




Fig. 3. — Application of the spectroscopic constraints on the hght ratio in B and V to 
determine k, the ratio of the radii of the components. The sohd curves are derived from 
hght curve solutions based on our CTIO photometry with k held fixed. The arrow with the 
error bar marks the resulting weighted average of k. 
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Fig. 4. — Light curves obtained at CTIO in the B and V bands along with the fitted solutions 
from EBOP. The lower panels expand the region near the primary and secondary minima 
for the V band. 
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Fig. 5. — Light curves in UBVR based on the observations by Zakirov (1997), along with 
our fits. 
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Fig. 6. — Comparison with stellar evolution models for both components in the logg^- 
logTeff plane. Evolutionary tracks for two different compositions are shown, giving good 
fits assuming GG Ori is on the main sequence (MS). The uncertainty in the placement of 
the tracks due to errors in the measured masses is indicated with the error bar in the lower 
left. The full extent of the main sequence band is shown in the inset, emphasizing the fact 
that the system is very close to the ZAMS. 
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Fig. 7. — Best fit evolutionary tracks assuming GG Ori is still in the prc-main sequence 
(PMS) phase. The uncertainty in the placement of the tracks due to errors in the measured 
masses is indicated with the error bar in the lower right. 
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Fig. 8. — Evolutionary track for a mass equal to that of the primary in GG Ori and a 
composition close to solar {Z = 0.020, Y = 0.2729). Also shown are all other binary systems 
that have accurate absolute dimensions available in which at least one component has a 
similar mass. The 8 stars represented differ in mass by less than 1%, and all are adequately 
fit by this model. 
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Fig. 9. — Critical values of the surface gravity at which the stars in GG Ori reach 
synchronization and circularization, as a function of orbital period. Curves for two different 
tidal mechanisms are shown: Zahn (1992) and Tassoul & Tassoul (1997). Evolution carries 
the stars upwards. Our measurements of logg for the system are indicated with crosses. 
Both theories predict circularization has not yet occurred (as observed), but they disagree 
as to the synchronization. 
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Table 1. KPNO radial velocity measurements of GG Ori, and residuals 
from the final spectroscopic orbit. 



HJD 


IWa 


RVb 


(O-C)a 


(O-C)b 




(2,400,000+) 


(km s~^) 


(km s~^) 


(kms~^) 


(kms~^) 


Orbital phase 


45684.8406 


+91.5 


-56.9 


+3.99 


+0.53 


0.394 


46070.8049 


-60.8 


+89.1 


-0.99 


-2.37 


0.087 


46071.8576 


-79.4 


+ 108.3 


+2.60 


-0.21 


0.398 


47778.9875 


+82.3 


-54.6 


-1.12 


-3.71 


0.906 


47779.9970 


+98.3 


-70.3 


-0.41 


-0.05 


0.199 


47781.9962 


-85.1 


+ 110.9 


-1.45 


-0.83 


0.507 


48150.9916 


+96.3 


-72.1 


-2.42 


+ 1.86 


0.236 


48152.0011 


+71.8 


-43.1 


-0.79 


-2.71 


0.027 


51243.7258 


-87.7 


+ 118.1 


+0.96 


+2.88 


0.167 


51246.7258 


+60.8 


-33.1 


-0.36 


+2.02 


0.168 


51247.6287 


+94.6 


-65.4 


+0.47 


+0.03 


0.162 
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Table 2. CfA radial velocity measurements of GG Ori, and residuals 
from the final spectroscopic orbit. 



H.JD 


RV,\''' 






\\j yj)B 




(2 400 000+1 


(km ) 


(km ) 


(km ) 


(km s""'^) 




50144.6930 


-9.01 


+39.00 


+0.79 


+0.53 


0.394 


50407.9143 


-73.46 


+97.71 


-2.16 


-2.37 


0.087 


50409.9760 


-10.06 


+36.68 


-1.84 


-0.21 


0.398 


50439.8695 


+69.42 


-47.59 


-2.98 


-3.71 


0.906 


50441.8148 


-90.88 


+119.08 


-0.57 


-0.05 


0.199 


50443.8593 


+35.21 


-9.50 


-2.04 


-0.83 


0.507 


50448.6911 


-79.78 


+110.30 


-0.13 


+ 1.86 


0.236 


50473.8339 


-22.15 


+48.85 


+0.71 


-2.71 


0.027 


50474.7629 


-93.01 


+126.09 


+ 1.37 


+2.88 


0.167 


50494.6588 


-94.59 


+125.22 


-0.21 


+2.02 


0.168 


50739.9868 


-93.01 


+123.32 


+ 1.46 


+0.03 


0.162 


50770.8719 


+97.91 


-66.73 


+0.27 


+2.43 


0.819 


50772.9141 


-89.25 


+118.60 


+0.44 


+0.10 


0.127 


50792.9273 


-95.67 


+117.64 


-2.32 


-4.53 


0.145 


50799.8040 


-93.11 


+125.32 


+0.10 


+3.29 


0.182 


50799.8119 


-91.45 


+122.44 


+ 1.60 


+0.56 


0.183 


50823.8998 


+95.85 


-68.85 


-2.18 


+0.71 


0.816 


50825.8988 


-88.51 


+117.77 


-2.11 


+2.56 


0.117 


50826.7972 


-73.96 


+103.43 


-0.53 


+ 1.22 


0.253 


50826.8636 


-68.20 


+99.13 


+ 1.21 


+0.94 


0.263 


50826.8738 


-68.02 


+100.13 


+0.75 


+2.58 


0.264 


50859.7244 


-86.99 


+113.14 


-1.50 


-1.16 


0.218 


50916.6656 


+99.47 


-73.36 


+0.48 


-2.85 


0.804 


51089.0199 


+100.83 


-72.47 


+ 1.35 


-1.46 


0.795 


51117.9532 


-95.72 


+124.68 


-1.33 


+ 1.46 


0.158 


51121.9141 


+100.86 


-68.75 


+2.02 


+ 1.62 


0.755 


51148.8862 


+98.99 


-69.36 


+ 1.71 


-0.56 


0.822 


51177.7818 


-91.89 


+124.19 


+ 1.64 


+ 1.84 


0.180 


51241.6644 


+100.05 


-71.21 


+ 1.74 


-1.38 


0.813 


OlZD / .DzOo 


+96. zz 


— DO. zz 


+ 0.04 


— O.OZ 




51268.6168 


+87.27 


-55.87 


+2.77 


+0.13 


0.877 


51292.6444 


+36.14 


-5.90 


+ 1.51 


+0.14 


0.500 


51472.9868 


+88.29 


-61.96 


-3.13 


+0.97 


0.695 


51474.0126 


+89.66 


-61.47 


-2.73 


+2.44 


0.850 


51475.9391 


-94.76 


+123.06 


-2.13 


+ 1.60 


0.140 


51533.8139 


+88.15 


-58.59 


+0.54 


+0.53 


0.868 


51539.8100 


+101.68 


-69.24 


+2.09 


+ 1.88 


0.772 


51566.8594 


+93.02 


-62.84 


+0.84 


+0.86 


0.851 


51569.7008 


-63.96 


+92.21 


-1.54 


+ 1.03 


0.279 


51592.8330 


+93.60 


-69.04 


-5.86 


+ 1.95 


0.767 


51619.6948 


+95.01 


-68.30 


-2.74 


+0.98 


0.818 


51621.7149 


-89.20 


+114.03 


-0.89 


-3.10 


0.123 



Includes corrections described in the text. 
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Table 3. Spectroscopic orbital solutions for GG Ori. 



Element KPNO CfA Combined 

(1) (2) (3) (4) 



Adjusted quantities 

(days)^. 6.6314948 6.6314948 6.6314936 ± 0.0000017 

7(kms-i) +14.13 ± 0.35 +14.30 ± 0.21 +14.31 + 0.21 

ARV (km s-^)'^ • • • • • • -0.06 + 0.39 

Ka (kms-^) 96.80 + 0.71 97.17 + 0.34 97.10 + 0.30 

i^B (kms-i) 96.67 + 0.51 97.53 + 0.33 97.28 + 0.28 

e 0.2264 + 0.0039 0.2201 + 0.0030 0.2218 + 0.0022 

iOA (deg) 122.7 + 1.0 123.34 + 0.67 122.76 + 0.39 

w (deg cycle-i) • • • • • • 0.00061 + 0.00025 

Min I (H,JD-2,400,000)^. . 49,717.66253 49,717.66253 49,717.66259 + 0.00020 

Derived quantities 

Pa (days) • • • • • • 6.6315048 ± 0.0000042 

U (yr) • • • • • • 10700 + 4500 

T (HJD-2,400,000)'=. 49,717.999 + 0.025 49,718.073 + 0.012 49,718.0448 + 0.0078 

OA sin« (IQfi km) 8.597 + 0.073 8.644 + 0.032 8.634 + 0.028 

Ob sini (10^ km) 8.586 + 0.054 8.676 + 0.031 8.650 + 0.026 

asini (R0) 24.69 + 0.13 24.885 + 0.065 24.833 + 0.056 

Ma sin^ i (Mq) 2.297 + 0.036 2.358 + 0.019 2.341 ± 0.016 

Mb sin^ I (Mq) 2.300 + 0.043 2.349 + 0.020 2.337 + 0.017 

q = Mb/Ma 1.001 + 0.010 0.9963 + 0.0051 0.9982 + 0.0044 

Other quantities pertaining to the fit 

A^obs 11 42 53 

A^minl ■■■ ■■■ 24 

-^min II • • • ■ ■ ■ 32 

Time span (yr) 'J 15.2 4.0 70.1 

CTA (kms-i) 1.93 1.95 1.83/1.93 

(TB (kms-i) 1.38 1.91 1.33/1.94 



^Sidereal period and epoch fixed from §3 for the KPNO and CfA solutions, and adjusted for the combined solution. 

''Velocity offset in the sense (CfA-KPNO). 

''Periastron passage closest to the epoch of primary minimum. 

'^Includes times of minimum for the combined solution. 
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Table 4. Eclipse timings and residuals for GG Ori. 



HJD 
(2,400,000+) 


Year 


Eclipse^ 


Type^ 


Epoch'' 


(O-C) 

(days) 
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2 


pg 
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-0.0017 


12 


41961.5670 


1973.763 


2 


pg 


-1170.0 


-0.0548 


11 


42074.3590 


1974.072 


2 


pg 


-1153.0 


+0.0019 


11 


42452.3650 


1975.107 


2 


pg 


-1096.0 


+0.0130 


11 


42469.2960 


1975.153 


1 


V 


-1093.0 


(-0.1399) 


13 


42469.3720 


1975.154 


1 


V 


-1093.0 


-0.0639 


14 


42777.3060 


1975.997 


2 


V 


-1047.0 


+0.0111 


15 


45403.3370 


1983.186 


2 


V 


-651.0 


-0.0276 


16 


48590.3086 


1991.912 


1 


pe 


-170.0 


+0.0006 


17 


48606.3690 


1991.956 


2 


V 


-168.0 


-0.0048 


18 


48623.4634 


1992.003 


1 


V 


-165.0 


-0.0021 


18 


48911.4212 


1992.791 


2 


pe 


-122.0 


-0.0011 


17 


48948.4090 


1992.892 


1 


pe 


-116.0 


+0.0001 


17 


49355.7327 


1994.007 


2 


pe 


-55.0 


+0.0006 


19 


49366.1968 


1994.036 


1 


pe 


-53.0 


+0.0036 


17 


49717.6624 


1994.998 


1 


pe 


0.0 


-0.0002 


17 



^1 = primary eclipse; 2 = secondary eclipse. 

^pg — photographic, v — visual, pe — photoelectric. 

'^Counted from primary minimum in column 4 of Table 3. 



References for Table 4. 

(I) Hoffmeister 1934; (2) Kordylewski 1951; (3) Soloviev 1945; (4) Nikulina 1958; (5) Zessewitsch 1945; 
(6) Szafraniec 1952; (7) Szafraniec 1955; (8) Tshuprina 1957; (9) Szafranicc 1956; (10) Szafraniec 1958; 

(II) Busch 1975; (12) Winiarski 1972; (13) Braune et al. 1977; (14) Diethelm 1975; (15) Diethelm 1976; 
(16) Diethelm 1983; (17) Lacy et al. 1995; (18) Diethelm 1992; (19) Lacy & Fox 1994. 
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Table 5. Absolute photometry of GG Ori and comparison stars. 



Star V B-V U-B iV„bs 



GG Ori (1988-1990) 10.359 ± 0.009 0.532 ± 0.004 0.323 ± 0.005 4 

GG Ori (1993-1995) 10.385 ± 0.007 0.490 ± 0.003 0.326 ± 0.005 3 

BD-1°1013 8.802 ± 0.003 0.232 ± 0.001 -0.224 ± 0.004 2 

BD-0°1081 8.960 ± 0.004 0.350 ± 0.004 0.027 ± 0.002 3 

BD-0°1044 7.341 0.523 -0.026 1 
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Table 6. Differential photometric observations of GG Ori from CTIO in the B 
band, in the sense (variable— comparison). 



HJD 




B 


HJD 




B 


HJD 




B 


HJD 




B 


(2,400,000+) 


(mag) 


(2,400,000+) 


(mag) 


(2,400,000+) 


(mag) 


(2,400,000+) 


(mag) 


49339 


6586 .... 


2 


536 


49342 


5826 


2 


038 


49355 


5552 


1 


884 


49717 


6308 


2 


333 


49339 


6677 .... 


2 


567 


49342 


5859 


2 


036 


49355 


5594 


1 


897 


49717 


6353 


2 


382 


49339 


6806 .... 


2 


497 


49343 


5543 


1 


867 


49355 


5625 


1 


903 


49717 


6391 


2 


408 


49339 


6869 .... 


2 


444 


49.34.3 


5590 


1 


886 


49355 


5658 


1 


900 


49717 


6427 


2 


439 


49339 


6919 . . . . 


2 


405 


49.34.3 


56.37 


1 


889 


49355 


5909 


1 


956 


49717 


6462 


2 


474 


49339 


6966 .... 


2 


360 


49344 


6894 


1 


859 


49355 


5953 


1 


938 


49717 


6498 


2 


492 


49339 


7145 .... 


2 


221 


A r\ O A A 

49 344 


6961 


1 


876 


49355 


5986 


1 


966 


49717 


6533 


2 


524 


49339 


7190 ... . 


2 


184 


49344 


7014 


1 


866 


49355 


6020 


1 


974 


49717 


6569 


2 


545 


49339 


7239 .... 


2 


162 


49348 


6986 


1 


870 


49355 


6052 


1 


991 


49717 


6613 


2 


543 


49339 


7283 .... 


2 


116 


49348 


7052 


1 


871 


49355 


6274 


2 


058 


49717 


6648 


2 


554 


49339 


7339 .... 


2 


099 


49348 


7102 


1 


873 


49355 


6318 


2 


059 


49717 


6684 


2 


543 


49339 


7385 .... 


2 


067 


49348 


8059 


1 


843 


49355 


6353 


2 


088 


49717 


6721 


2 


527 


49339 


7446 .... 


2 


023 


49348 


8101 


1 


825 


49355 


6388 


2 


110 


49717 


6757 


2 


486 


49339 


7487 .... 


2 


010 




5506 


1 


9.31 


49355 


6438 


2 


127 


49717 


6793 


2 


464 


49339 


7527 .... 


1 


990 


49349 


5606 


1 


887 


49355 


6487 


2 


140 


49717 


6840 


2 


432 


49339 


7607 .... 


1 


941 


49349 


5673 


1 


868 


49355 


6730 


2 


260 


49717 


6872 


2 


402 


49339 


7684 .... 


1 


923 


49349 


6409 


1 


881 


49355 


6775 


2 


274 


49717 


6907 


2 


371 


49339 


7754 .... 


1 


894 


49349 


6454 


1 


889 


49355 


6807 


2 


304 


49717 


6943 


2 


341 


49339 


7820 .... 


1 


882 


49349 


6488 


1 


886 


49355 


6840 


2 


316 


49717 


6978 


2 


296 


49339 


7886 .... 


1 


868 


49349 


6533 


1 


882 


49355 


6874 


2 


343 


49717 


7037 


2 


262 


49339 


7944 .... 


1 


858 


49349 


6566 


1 


889 


49355 


6907 


2 


357 


49717 


7073 


2 


223 


49339 


8041 .... 


1 


852 


49349 


6610 


1 


878 


49355 


6950 


2 


389 


49717 


711 1 


2 


223 


49339 


8086 .... 


1 


855 


49349 


6642 


1 


876 


49355 


G9S3 


2 


386 


49717 


7146 


2 


179 


49339 


8138 . . . . 


1 


854 


49.349 


6683 


1 


861 


49355 


7017 


2 


419 


49717 


7192 


2 


149 


49339 


8181 .... 


1 


854 


49349 


6716 


1 


866 


49355 


7050 


2 


417 


49717 


7228 


2 


136 


49339 


8225 .... 


1 


856 


49349 


7644 


1 


862 


49355 


7086 


2 


454 


49717 


7264 


2 


104 


49339 


8267 .... 


1 


837 


49349 


7688 


1 


862 


49355 


7132 


2 


467 


49717 


7304 


2 


071 


49339 


8331 .... 


1 


881 


49349 


7722 


1 


864 


49355 


7166 


2 


465 


49717 


7340 


2 


059 


49340 


5.368 .... 


1 


829 


49349 


7755 


1 


863 


49355 


7201 


2 


493 


49717 


7376 


2 


028 


49340 


5453 . . . . 


1 


926 


49349 


7812 


1 


857 


49355 


7234 


2 


505 


49717 


7412 


2 


019 


49340 


5520 .... 


1 


908 


49349 


7847 


1 


856 


49355 


7270 


2 


509 


49717 


7447 


1 


992 


49340 


5578 .... 


1 


919 


49.350 


58G7 


1 


897 


49355 


7305 


2 


506 


49717 


7485 


1 


973 


49340 


5639 .... 


1 


896 


49350 


5902 


1 


879 


49355 


7341 


2 


514 


49717 


7563 


1 


940 


49340 


5680 .... 


1 


884 


49350 


5933 


1 


892 


49355 


7376 


2 


502 


49717 


7626 


1 


915 


49340 


5736 .... 


1 


889 


49350 


6623 


1 


891 


49355 


7421 


2 


502 


49717 


7681 


1 


899 


49340 


5776 .... 


1 


884 


49350 


6686 


1 


894 


49355 


7455 


2 


488 


49717 


7750 


1 


866 


49340 


5821 .... 


1 


894 


49350 


6762 


1 


880 


49355 


7494 


2 


477 


49717 


7799 


1 


863 


49340 


5862 .... 


1 


889 


49.350 


7574 


1 


861 


49355 


7536 


2 


453 


49717 


7862 


1 


859 


49340 


5915 . . . . 


1 


897 


49350 


7620 


1 


855 


49355 


7571 


2 


436 


49717 


7972 


1 


849 


49340 


5955 .... 


1 


878 


49350 


7654 


1 


859 


49716 


5883 


1 


878 


49720 


5503 


2 


058 


49340 


6794 .... 


1 


870 


49350 


8027 


1 


861 


49716 


5977 


1 


884 


49720 


5547 


2 


154 


49340 


6857 


1 


869 


49350 


8074 


1 


838 


49716 


6064 


1 


885 


49720 


5556 


2 


138 


49340 


6898 


1 


869 


49350 


8114 


1 


837 


49716 


7719 


1 


873 


49720 


5601 


2 


130 


49340 


6943 


1 


868 


49352 


7154 


1 


867 


49716 


7784 


1 


847 


49720 


5646 


2 


100 


49340 


6982 


1 


878 


49352 


7214 


1 


866 


49716 


7831 


1 


850 


49720 


5680 


2 


097 


49340 


7018 


1 


878 


49352 


7266 


1 


858 


49717 


5705 


1 


971 


49720 


5713 


2 


053 


49340 


7573 


1 


868 


49352 


7598 


1 


846 


49717 


5757 


1 


991 


49720 


5755 


2 


047 


49340 


7661 


1 


857 


49352 


7644 


1 


861 


49717 


5799 


2 


010 


49720 


5788 


2 


036 


49340 


7726 


1 


866 


49352 


7684 


1 


860 


49717 


5846 


2 


029 


49720 


5823 


2 


017 


49341 


7584 


1 


873 


49352 


7721 


1 


859 


49717 


5883 


2 


048 


49720 


5859 


2 


000 


49341 


7664 


1 


860 


49353 


5467 


1 


883 


49717 


5920 


2 


066 


49720 


5891 


2 


000 


49341 


7731 


1 


860 


49353 


5515 


1 


899 


49717 


5958 


2 


085 


49720 


5934 


1 


984 


49342 


5482 


2 


194 


49353 


5549 


1 


890 


49717 


6005 


2 


127 


49720 


5969 


1 


967 


49342 


5544 


2 


155 


49353 


5585 


1 


892 


49717 


6041 


2 


151 


49720 


6021 


1 


956 


49342 


5590 


2 


158 


49354 


6393 


1 


874 


49717 


6080 


2 


160 


49720 


6071 


1 


949 


49342 


5635 


2 


116 


49354 


6439 


1 


869 


49717 


6116 


2 


192 


49720 


7236 


1 


880 


49342 


5668 


2 


114 


49354 


6474 


1 


868 


49717 


6164 


2 


243 


49720 


7289 


1 


874 


49342 


5714 


2 


096 


49355 


5410 


1 


891 


49717 


6198 


2 


253 


49720 


7326 


1 


867 


49342 


5747 


2 


071 


49355 


5468 


1 


884 


49717 


6233 


2 


280 


49720 


7814 


1 


865 


49342 


5791 


2 


062 


49355 


5521 


1 


890 


49717 


6271 


2 


320 


49720 


7873 


1 


858 


49720 


7913 


1 


854 


49722 


6534 


1 


882 


49722 


7456 


1 


855 


49727 


5643 


1 


906 


49721 


7524 


1 


849 


49722 


6587 


1 


870 


49726 


6865 


1 


878 


49727 


5700 


1 


908 


49721 


7590 


1 


857 


49722 


7343 


1 


874 


49726 


6927 


1 


891 










49721 


7643 


1 


849 


49722 


7378 


1 


854 


49726 


6981 


1 


883 










49722 


6461 


1 


883 


49722 


7423 


1 


852 


49727 


5591 


1 


907 
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Table 7. Differential photometric observations of GG Ori from CTIO in the V 
band, in the sense (variable— comparison). 



HJD 




V 


HJD 




V 


HJD 




V 


HJD 




V 


(2,400,000+) 


(mag) 


(2,400,000+) 


(mag) 


(2,400,000+) 


(mag) 


(2,400,000+) 


(mag) 


49339 


6590 .... 


2 


236 


49342 


5830 


1 


754 


49355 


5555 


1 


605 


49717 


6312 


2 


074 


49339 


6681 .... 


2 


291 


49342 


5863 


1 


739 


49355 


5597 


1 


607 


49717 


6357 


2 


096 


49339 


6810 . . . . 


2 


224 


49343 


5546 


1 


579 


49355 


5629 


1 


626 


49717 


6395 


2 


144 


49339 


6873 .... 


2 


164 


49.34.3 


5594 


1 


606 


49355 


5662 


1 


620 


49717 


6430 


2 


161 


49339 


6923 .... 


2 


124 


49.34.3 


5641 


1 


609 


49355 


5913 


1 


672 


49717 


6465 


2 


197 


49339 


6970 .... 


2 


088 


49344 


6898 


1 


604 


49355 


5957 


1 


683 


49717 


6502 


2 


220 


49339 


7149 .... 


1 


947 


A r\ O A A 

49344 


6965 


1 


593 


4935.5 


5990 


1 


68S 


49717 


6537 


2 


242 


49339 


7194 ... . 


1 


916 


49344 


7018 


1 


600 


49355 


6024 


1 


691 


49717 


6573 


2 


268 


49339 


7243 .... 


1 


878 


49348 


6990 


1 


601 


49355 


6055 


1 


709 


49717 


6617 


2 


277 


49339 


7287 .... 


1 


856 


49348 


7056 


1 


589 


49355 


6278 


1 


790 


49717 


6652 


2 


277 


49339 


7343 .... 


1 


823 


49348 


7106 


1 


585 


49355 


6322 


1 


808 


49717 


6689 


2 


267 


49339 


7389 .... 


1 


792 


49348 


8063 


1 


586 


49355 


6357 


1 


813 


49717 


6725 


2 


235 


49339 


7450 .... 


1 


754 


49.348 


8105 


1 


572 


49355 


6392 


1 


826 


49717 


6761 


2 


211 


49339 


7491 .... 


1 


741 


A r\o A Ci 

4y.34y 


5509 


1 


632 


49355 


6442 


1 


841 


49717 


6797 


2 


181 


49339 


7531 .... 


1 


709 


49349 


.5610 


1 


602 


49355 


6491 


1 


874 


49717 


6844 


2 


143 


49339 


7611 .... 


1 


679 


49349 


.5677 


1 


601 


49355 


6734 


1 


986 


49717 


6876 


2 


124 


49339 


7688 .... 


1 


648 


49349 


6413 


1 


591 


49355 


6779 


2 


016 


49717 


6911 


2 


088 


49339 


7758 .... 


1 


626 


49349 


6458 


1 


596 


49355 


6811 


2 


025 


49717 


6947 


2 


059 


49339 


7824 .... 


1 


601 


49349 


6492 


1 


609 


49355 


6844 


2 


042 


49717 


6982 


2 


027 


49339 


7890 .... 


1 


595 


49349 


6537 


1 


600 


49355 


6878 


2 


061 


49717 


7041 


1 


983 


49339 


7948 .... 


1 


599 


49349 


6570 


1 


605 


49355 


6911 


2 


074 


49717 


7077 


1 


957 


49339 


8045 .... 


1 


.591 


49349 


6614 


1 


585 


49355 


69.54 


2 


097 


49717 


7115 


1 


937 


49339 


8090 .... 


1 


594 


49349 


6646 


1 


596 


49355 


G9S7 


2 


122 


49717 


7150 


1 


913 


49339 


8142 .... 


1 


596 


A r\o A (\ 

4y.34y 


G687 


1 


591 


49355 


7020 


2 


137 


49717 


7196 


1 


881 


49339 


818.5 . . . . 


1 


594 


40349 


6720 


1 


595 


49355 


7054 


2 


164 


49717 


7232 


1 


855 


49339 


8229 .... 


1 


582 


49349 


7648 


1 


598 


49355 


7090 


2 


175 


49717 


7268 


1 


832 


49339 


8271 .... 


1 


590 


49349 


7692 


1 


586 


49355 


7136 


2 


193 


49717 


7308 


1 


805 


49339 


8335 .... 


1 


650 


49349 


7726 


1 


589 


49355 


7170 


2 


219 


49717 


7344 


1 


776 


49340 


5.372 .... 


1 


577 


49349 


7759 


1 


580 


49.355 


7205 


2 


221 


49717 


7380 


1 


766 


49340 


5456 .... 


1 


621 


49349 


7816 


1 


583 


4935.5 


7237 


2 


243 


49717 


7416 


1 


736 


49340 


5523 .... 


1 


609 


49349 


7850 


1 


582 


4935.5 


7274 


2 


238 


49717 


7451 


1 


724 


49340 


5582 .... 


1 


610 


49350 


5871 


1 


604 


49355 


7309 


2 


2.37 


49717 


7489 


1 


707 


49340 


5643 .... 


1 


618 


49350 


5906 


1 


602 


49355 


7345 


2 


245 


49717 


7567 


1 


680 


49340 


5684 .... 


1 


611 


49350 


5937 


1 


614 


49355 


7380 


2 


2.36 


49717 


7630 


1 


651 


49340 


5739 .... 


1 


611 


49350 


6627 


1 


578 


49355 


7425 


2 


232 


49717 


7685 


1 


627 


49340 


5780 .... 


1 


607 


49350 


6690 


1 


627 


49355 


7459 


2 


229 


49717 


7754 


1 


607 


49340 


5825 .... 


1 


611 


49350 


6766 


1 


600 


49355 


7498 


2 


213 


49717 


7803 


1 


591 


49340 


5866 .... 


1 


605 


49.350 


7578 


1 


603 


49355 


7539 


2 


194 


49717 


7866 


1 


595 


49340 


5919 . . . . 


1 


595 


49350 


7624 


1 


586 


49355 


7574 


2 


174 


49717 


7976 


1 


.585 


49340 


5959 .... 


1 


606 


49350 


7658 


1 


600 


49716 


.5887 


1 


602 


49720 


5505 


1 


819 


49340 


6798 .... 


1 


600 


49350 


8031 


1 


579 


49716 


.5981 


1 


598 


49720 


5548 


1 


844 


49340 


6861 


1 


596 


49350 


8078 


1 


559 


49716 


6068 


1 


597 


49720 


5560 


1 


860 


49340 


6902 


1 


595 


49350 


8118 


1 


586 


49716 


7723 


1 


597 


49720 


5606 


1 


825 


49340 


6947 


1 


592 


49352 


7158 


1 


590 


49716 


7788 


1 


590 


49720 


5650 


1 


814 


49340 


6986 


1 


599 


49352 


7218 


1 


599 


49716 


7835 


1 


585 


49720 


5684 


1 


791 


49340 


7022 


1 


600 


49352 


7269 


1 


587 


49717 


5709 


1 


691 


49720 


5717 


1 


780 


49340 


7577 


1 


598 


49352 


7602 


1 


595 


49717 


5761 


1 


716 


49720 


5759 


1 


771 


49340 


7665 


1 


590 


49352 


7648 


1 


595 


49717 


5802 


1 


747 


49720 


5792 


1 


756 


49340 


7730 


1 


596 


49352 


7689 


1 


591 


49717 


5850 


1 


753 


49720 


5827 


1 


753 


49341 


7588 


1 


593 


49352 


7725 


1 


585 


49717 


5887 


1 


778 


49720 


5863 


1 


734 


49341 


7668 


1 


591 


49353 


5471 


1 


604 


49717 


5924 


1 


802 


49720 


5895 


1 


717 


49341 


7735 


1 


591 


49353 


5519 


1 


608 


49717 


5962 


1 


815 


49720 


5938 


1 


706 


49342 


5486 


1 


917 


49353 


5554 


1 


612 


49717 


6009 


1 


849 


49720 


5973 


1 


698 


49342 


5548 


1 


875 


49353 


5589 


1 


606 


49717 


6045 


1 


865 


49720 


6025 


1 


677 


49342 


5594 


1 


859 


49354 


6397 


1 


602 


49717 


6084 


1 


897 


49720 


6075 


1 


661 


49342 


5638 


1 


829 


49354 


6443 


1 


594 


49717 


6120 


1 


925 


49720 


7240 


1 


609 


49342 


5672 


1 


838 


49354 


6478 


1 


602 


49717 


6168 


1 


960 


49720 


7293 


1 


599 


49342 


5718 


1 


802 


49355 


5414 


1 


610 


49717 


6201 


1 


983 


49720 


7329 


1 


588 


49342 


5751 


1 


787 


49355 


5472 


1 


602 


49717 


6237 


2 


003 


49720 


7818 


1 


591 


49342 


5795 


1 


774 


49355 


5525 


1 


604 


49717 


6275 


2 


038 


49720 


7877 


1 


590 


49720 


7917 


1 


578 


49722 


6538 


1 


595 


49722 


7460 


1 


584 


49727 


5647 


1 


624 


49721 


7528 


1 


591 


49722 


6591 


1 


593 


49726 


6869 


1 


602 


49727 


5704 


1 


616 


49721 


7594 


1 


595 


49722 


7347 


1 


594 


49726 


6931 


1 


605 










49721 


7647 


1 


592 


49722 


7382 


1 


583 


49726 


6985 


1 


609 










49722 


6465 


1 


592 


49722 


7427 


1 


585 


49727 


5594 


1 


622 
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Table 8. Light curve solutions for GG Ori based on our CTIO observations. 



General Fits Adopted Fits 



Parameter B V B V 

(1) (2) (3) (4) (5) 



Jb 0.991 ± 0.030 0.992 ± 0.018 1.002 ± 0.008 1.002 ± 0.005 

rA 0.074 ± 0.003 0.073 ± 0.002 0.0746 ± 0.0003 

k = rB/rA^- 1.02 ± 0.08 1.04 ± 0.05 0.988 

« (deg) 89.20 ± 0.14 89.31 ± 0.11 89.24 ± 0.08 

ecoscj -0.1187 ± 0.0003 -0.1193 ± 0.0002 -0.1191 ± 0.0003 

esinw +0.210 ± 0.009 +0.198 ± 0.005 +0.200 + 0.007 

u = ua = ub .... 0.59 + 0.12 0.48 + 0.08 0.55 + 0.05 0.55 + 0.03 

La 0.492 + 0.043 0.482 + 0.028 0.504 + 0.006 0.504 + 0.005 

(J (mag) 0.01414 0.00879 0.01408 0.00884 

A^obs 256 255 256 255 



^In the final fits k was held fixed to the value derived from spectroscopic constraints described in the text. 
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Table 9. Light curve solutions based on the UBVR observations by Zakirov (1997). 



Parameter UBVR 



Restricted fits 

Jb 0.98 ± 0.04 1.006 ± 0.012 0.991 ± 0.008 1.006 ± 0.012 

ta 0.076 ± 0.002 0.0734 ± 0.0007 0.0746 ± 0.0005 0.0754 ± 0.0007 

k = rB/rA''- 0.988 0.988 0.988 0.988 

i (deg) 89.5 ± 0.3 89.43 ± 0.09 89.35 ± 0.07 89.24 ± 0.09 

ecoscj -0.1197 ± 0.0013 -0.1202 ± 0.0002 -0.1197 ± 0.0002 -0.1194 ± 0.0002 

esinw +0.180 ± 0.024 +0.179 ± 0.008 +0.178 + 0.006 +0.196 + 0.008 

u = UA = UB.... 0.49 (fixed) 0.36 + 0.11 0.41 + 0.09 0.33 + 0.13 

La 0.510 + 0.034 0.503 + 0.011 0.507 + 0.008 0.503 + 0.011 

cr (mag) 0.02499 0.01837 0.01365 0.01774 

A^obs 175 240 242 238 

Adopted fits 

Jb 0.964 + 0.016 1.004 + 0.009 0.990 + 0.007 1.009 + 0.009 

ta 0.0746 + 0.0005 

k = rB/rA''. 0.988 

i (deg) 89.35 + 0.07 

ecoscj -0.1198 + 0.0004 

esinw +0.184 + 0.008 

u = UA = UB.... 0.49 (fixed) 0.41 + 0.07 0.41 + 0.05 0.23 + 0.07 

La 0.513 + 0.011 0.504 + 0.008 0.507 + 0.008 0.502 + 0.008 

a (mag) 0.02519 0.01848 0.01361 0.01781 

iVobs 175 240 242 238 



^Held fixed in all solutions, based on the spectroscopic constraints described in the text. 
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Table 10. Weighted average light elements for GG Ori. 



Parameter U B V R 



Jb 0.964 ± 0.016 1.003 ± 0.006 0.998 ± 0.004 1.009 ± 0.009 

rA 0.0746 ± 0.0010 

rs^ 0.0737 ± 0.0010 

r = r^+rB^.... 0.1483 ± 0.0019 

k = rB/rA 0.988 ± 0.009 

i (deg) 89.30 ± 0.10 

ecoscj -0.1194 ± 0.0002 

esinw +0.193 ± 0.005 

u = UA = UB .... 0.49 (fixed) 0.50 ± 0.04 0.51 ± 0.03 0.23 ± 0.07 

La 0.513 ± 0.011 0.504 ± 0.005 0.505 ± 0.004 0.502 ± 0.008 



'^See text for error estimates. 
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Table 11. Physical parameters of GG Ori. 



Parameter 



Mass (Mq) 

Radius (Rq) 

logs 

p (gr cm~^) 

vsini (kms^i) (KPNO) . . . . 

vsini (kms^^) (CfA) 

^sync (kms~^) (circular) . . . . 
t^sync (kms^^) (eccentric)... 

« (R©) 

logTeff 

logL/L© 

Mboi (mag) 



My (mag) . . . 

Lb /La 

m — M (mag) 
Dist (pc) .... 



Primary Secondary 



2.342 


± 


0.016 


2.338 


± 


0.017 


1.852 


± 


0.025 


1.830 


± 


0.025 


4.272 


± 


0.012 


4.282 


± 


0.012 


0.519 


± 


0.021 


0.537 


± 


0.022 


24 


± 


2 


23 


± 


2 


16 


± 


1 


16 


± 


1 


14.1 


± 


0.2 


14.0 


± 


0.2 


22.8 


± 


0.3 


22.5 


± 


0.3 






24.835 


± 0.056 






3.9978 


± 


0.0088 


3.9978 


± 


0.0088 


1.480 


± 


0.037 


1.470 


± 


0.037 


1.020 


± 


0.093 


1.046 


± 


0.092 


3.974 


± 


0.003 


3.974 


± 


0.003 


1.176 


± 


0.038 


1.202 


± 


0.039 






0.98 


± 0.11 










8.20 


± 0.10 










438 


± 20 
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Table 12. Evolutionary age of the components of GG Ori in 
the main-sequence (MS) and pre-main sequence (PMS) scenarios, 

in Myr. 





Chemical composition 


Primary 


Secondary 


Average 


z = 


0.020, F - 0.2729 (MS) 


93 ±27 


72 ±27 


82± 19 


z = 


0.021, y = 0.2810 (MS) 


75 ± 29 


52 ±32 


65 ±21 


z = 


0.022, F = 0.2860 (PMS) . . . . 


7.26 


8.23 


7.7 



